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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 206k 

EFFECT OF ASPECT RATIO ON THE AIR FORCES AND MOMENTS OF 
HARMONICALLY OSCILLATING THIN RECTANGULAR 
WINGS IN SUPERSONIC POTENTIAL FLOW 
By Charles E. Watkins 

SUMMARY 


This paper treats the effect of aspect ratio on the air forces and 
moments of an oscillating flat rectangular wing in supersonic potential 
flow. The linearized velocity potential for the wing undergoing 
sinusoidal torsional oscillations simultaneously with sinusoidal 
vertical translations is derived in the form of a power series in terms 
of a frequency parameter. The series development is such that the 
differential equation for the velocity potential is satisfied to the 
required power of the frequency parameter considered and the linear 
boundary conditions are satisfied exactly. The method of solution can 
be utilized for other plan forms, that is plan forms for which certain 
steady- state solutions are known. 

Simple, closed expressions that include the reduced frequency to - 
the third power, which is thought to be sufficient for most practical 
applications, are given for the velocity potential, the components of 
total force and moment coefficients, and the components of chordwise 
section force and moment coefficients . The components of total force 
and moment coefficients indicate the over -all effect of aspect ratio 
on these quantities 3 however, the components of chordwise coefficients 
yield more information because they account for the spanwise distribu- 
tion of aerodynamic loading of a rectangular wing and may therefore be 
more useful for flutter calculations. It is found that the components 
of force and moment coefficients for a small-aspect-ratio wing may 
deviate considerably from those for an infinite -aspect-ratio wing. 
Thickness effects which may alter some of the conclusions are not taken 
into account in the analysis. Results of some selected calculations ar e 
presented in several figures and discussed. 
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INTRODUCTION 


The. effect of aspect ratio on the single-degree torsional 
instability of a finite rectangular wing oscillating in a supersonic 
stream was treated in reference 1 by expanding, in powers of the fre- 
quency of oscillation, the linearized velocity potential developed in 
reference 2. Since only slow oscillations were considered pertinent to 
single -degree torsional instability, terms in the expansion involving 
the frequency of oscillation to powers higher than the first were not 
considered. 

In the present paper the expanded linearized velocity potential is 
used to study the effect of aspect ratio on the air forces and moments 
of an oscillating, thin, flat, finite, rectangular wing when higher 
powers of the frequency of oscillation are taken into account. The 
motions considered are sinusoidal torsional oscillations about a span- 
wise axis taken simultaneously with sinusoidal vertical translations of 
this axis . The velocity potential is developed by superpositions of 
sources and doublets, so as to include all powers of the frequency of 
oscillations up to any desired power. Simple, closed expressions are 
given for the velocity potential, components of the total force and 
moment coefficients, and components of the chordwise section force and 
moment coefficients involving powers of the frequency up to and 
including the third power. Extension of the results to include higher 
powers of the frequency is straightforward. 

A recent publication, reference 3 , that became available after this 
investigation was completed, is partly devoted to the treatment of a 
rectangular wing undergoing the same types of harmonic motions as those 
considered herein. The velocity potential is determined in the form of 
a double integral, by application of the Fourier transform to the 
boundary-value problem for this potential, and expressions for forces 
and moments are given in terms of this double integral. The reduction 
of the integral expressions of reference 3 to forms desirable for 
flutter calculations, that is, chordwise section forces and moments, 
is not given. 


SYMBOLS 


0 disturbance -velocity potential 

rectangular coordinates attached to wing moving in 
negative x-direction 


x,y,z 
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M 


2m 


w(x, yi ,t) 


x 0 


rectangular coordinates used, to represent space location 
of sources or doublets in xy-plane 

function defining mean ordinates of any chordvrise section 
of wing such as y = y^ as shown in figure 1 

vertical velocity at surface of wing along chordwise 
section at y = y^ 

abscissa of axis of rotation of wing (elastic axis) as 
shown in figure 1 - 


t 


time 


h 

*0 

a 

a 0 

h,a 

V 


vertical displacement of axis of rotation 

amplitude of vertical displacement of axis of rotation, 
positive downward 

angle of attack 

amplitude of angular displacement about axis of rotation, 
positive leading edge up 

time derivative of h and a, respectively 

velocity of main stream 


c 

M 


P = \/m 2 - 1 

t 1j t 2^i^2 

w(i,ti) 

w(t) 

CO 


_ M 2 o> 
“ = — 

Vp^ 


velocity of sound 

free -stream Mach number (V/c) 

functions defined with equation (7) 

function used to represent space variation of source and 
doublet strengths 

function used to represent time variation of source and 
doublet strengths 

frequency of oscillation 


reduced frequency (ccib/V) 
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R = 3^ 

a nm 


f l' f 2> f 3 


D i 


Gn 

P 

Ap 


b 


ni)0i2 - n) 

represents functions of 55, x ; and M, defined in 
equation ( 15 ) 

represent functions of x ; Xq, and <5, defined in 
equation ( 19 ) 

function used to denote doublet distributions (see 
equation (21) ) 

function defined in equation ( 28 ) 

function defined in equation ( 29 ) 

density- 

local pressure difference measured positive downward,, 
defined in equation ( 31 ) 

half -chord 


s 


half -span 


A 

P 

L 1> L 2-» L 3-> L 4 


Ml, Mg, M3, Mi,. 


P 


aspect ratio (s/b) 

total force acting on wing' defined in equation ( 32 ) 

components of total force coefficients, defined in 
equation ( 35 ) 

total moment acting on wing, defined in equation ( 36 ) 

components of total moment coefficients, defined in 
equation ( 38 ) 

section force (total force at any spanwise station), 
defined in equation (39) 


1*1 ,1*2,13 , 1*4 


components of section force coefficients, defined in 
equations (4l) and (42) 


Ma 


section moment (total moment at any spanwise station), 
defined in equation (40) 


Mi,M2,M 3 ,M4 


F n> G n 


components of section moment coefficients, defined in 
equations (43) and .(44) 

functions related to F n and G n , defined in appendix 
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ANALYSIS 

Boundary-Value Problems for Velocity Potentials 


Consider a thin flat rectangular wing moving at a constant super- 
sonic speed in a chordwise direction normal to its leading edge as 
shown in figure 1. The boundary -value problems for the velocity 
potential for such a wing may be conveniently classified' into two types 
associated with the nature of the flow over different portions of the 
wing. On the portion of the wing between the Mach cones emanating from 
the foremost point of each tip (region IT in fig. 1(a)) there is no 
interaction between the flow on the upper and lower surfaces of the 
wing. The type of boundary -value problem for this portion of the wing - 
is referred to herein as "purely supersonic" and the velocity potential 
for region IT is denoted by On portions of the wing within the 

Mach cones emanating from the foremost point of each tip (regions Tq, 
T 2 , and T 3 in fig. 1 (a) ), there is interaction between the flow on 
the upper and lower surfaces of the wing. The type of boundary-value 
problem for these portions of the wing is referred to as "mixed 
supersonic" and the velocity potentials for these regions are designated 
by 0qq, fa 2 ) an d 0 T 3 > respectively. The complete velocity potential 


at a point may then be expressed as 


0M 0T-]/ 0T 2 > 


or 



according • 


to the region that contains the point. 


As customary in linear theory, as applied to thin flat surfaces, 
the boundary conditions are to be ultimately satisfied by the velocity 
potentials at the projection of the wing onto a plane (the xy-plane) 
with respect to which all deflections are considered small and which 
lies parallel to the free-stream direction. Thickness effects - are not 
taken into account} hence, the velocity potentials are associated only 
with’ conditions that yield lift and are consequently antisymmetrical 
with respect to the plane of the projected wing. It is therefore 
necessary to consider the potentials at only one surface, upper or 
lower, of the projected wing. The upper surface is chosen for this 
analysis. 


The differential equation for the propagation of small disturbances 
that must be. satisfied by the velocity potentials is (when referred to 
a rectangular coordinate system x,y, z with the xy-plane coincident 
with the reference plane and moving uniformly in the negative 
x-direction, fig. 1 ) 
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The houndary conditions that must he satisfied hy the velocity potential 
are: (a) In regions T^_, Tg, T^, and N the flow must he tangent to the 

surface of the wing or 



wCx^y-j^t) = V 


chT + 


dt 


( 2 ) 


where is the vertical displacement of the ordinates of the surface 

of any chordwise section of the wing (see fig. 1(h)). (h) In regions Tq 

and T2 the pressure must fall to zero along the wing tips and remain 
zero in the portion of the Mach cones emanating from the foremost 
po int s of the wing tips not occupied hy the wing.' (Another condition, 
that the potential must he zero ahead of the wing and in the regions 
off the wing adjacent to the Mach cones emanating from the foremost 
points of the tips, is automatically satisfied hy the type of source 
and doublet synthesis employed in the solutions.) 

For the particular case of a wing independently performing small 
sinusoidal torsional oscillations of amplitude |oq| and frequency ca 

about some spanwise axis Xq and small sinusoidal vertical translations 
of amplitude |1 iq| and frequency (a, the equation of is 

= e^j^x - xq) + hq] = a(x - xq) + h (3) 

Substituting this expression for 2^ into equation (2) gives 

w(x,yp,t) = Va + a(x - xq) + h (4) 

The velocity potential may thus he expressed as the sum of separate 
effects due to position and motion of the wing associated with the 
individual terms in equation (4) as 

0 = + 0d + K ' (5) 


Derivation of 0^ 

The boundary -value problem in the purely supersonic region (fig. 2(a)) 
is the same as that for the two-dimensional wing treated in reference 4. 
This problem is there shown to he satisfied hy a distribution of sources 
referred to, in this case, as moving sources because of the uniform 
motion} that is. 
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0n(x,y,z,t) 



w (^Tl)01 An cli 


( 6 ) 


In equation (6) W(l,q) represents the space variation of source 
strength and must be evaluated in accordance with the individual te rns 
of equation ( 4 ) and is 'the potential of a moving source situated 

at the point (i,T|,0) that may be expressed as 


v(t - t x ) + w (t - Tg) 

- Hi)(il 2 - n) 


( 7 ) 


where w(t) is the time variation of source strength and the symbols 
with subscripts appearing in equation (7) are defined as 


M(x - I) _ \|(t) - m)(ng - 4) 


_ M(x - |) \/(n - ^i) (-ng - 4) 

2 CP 2 x P c 


1 /7 172 n 2 2 

Tii “ y - p\r x " ^ - 0 z 


42 = y 


■\l( x - 6 > 


2 - P 2 z 2 


The time variation of source strength w(t) for harmonic oscillations 
may be written as 


io>t 


w(t) = e 

The numerator in equation (7) thus becomes 

w(t - T X ) + w(t - Tg) = gi a> ( t-T l) + e i0)(t-Tg) 


(8) 


icut -io)-~ — t 2 “ T l , 

= 2e e 2 cos a> 5 (9) 


o 
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Substituting equations (7) and ( 9 ) into equation (6) yields 
Mw , t) - . w(M)e ' 15U ' e ’ -(!*)*> *1 (10) 
where, for briefness, 

_ ccM 

05 “ 2 - 2 

cr VP 

and 

R =\/(x - I) 2 - P 2 (y - T]) 2 - p 2 z 2 = 3^(n - t1 1 )(ti 2 - ^ 

The values of W( £,t| ) associated with the different terms of equa- 
tion (4) are, in the order in which they are used:- 


For h 


9 

^ iVp 2 co 

- ^ ho 

(11) 

For Va 

¥d,T]) = Voo 

(12) 

For a(x - xq) 

w( M) = P 2 ao( ^ - X 0 ) 
Mr 

(13) 

If any of the values 

of W(5 ,tj) given in equations' (ll) , 

(12),* and ( 13 ) 


is put into equation (10 ), the integration with respect to rj can be 
readily performed and the remaining integral evaluated as a series of 
Bessel functions. (See, for example, reference 4.) However, in order 
to be consistent with and to lead naturally to a succeeding part of the 
analysis the integrand is expanded into a Maclaurin's series with 
respect to o3. The expansion yields 
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icDt pX-| 3z p T]p FT 

^f(x,y,z,t) = - — J ^ J W(|,ri) [i 


a nl 


a 01 r + a ll r + • • • + 

~r~ + * • ' J + ( a 02^ + a i2^ + * • • + a^pl^R + . . . ) + . . . + 


3 2m-3 


3 2m-3 


, a O m R " + + * 


+ a 


v 1 ”® 2 "' 3 + • • •) + • • • 


<3.t| a| (14) 


where the coefficients a -n TO are f met ions of cd, x , and M; those 
coefficients involving cd, up to and including the third power, are: 


f£> 

O 

H 

II 

n . — UJ 

1 - 10XX ? 

c 

a ll = 

•- , -2 
1CD + CD X - - 

®21 = 

-2 3 

CD + 1CD J . 

2 2 ■ 

a 3l = 

—3 

ICD-* 

6 

II 

OJ 

o 

as 

_2 3 

_ 00 + 

2M 2 2M 2 

a 12 = 

--3 

XCD J 

2 

2&T 


_2 3 

2 lax 3 

+ ~r 


-3 

l /iw 


\ 


> 


X 


(15) 


J 

Observe the following identity- that is valid regardless of the highest 
power of cd considered and that will be of subsequent use, namely 


a 01 + xa "i 1 + • 


11 


+ x I1 a nl = 1 


(16) 
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It -will be noted in equation (l4) that the potential of a moving source 
when expanded in terms of the frequency appears as a series of terms 
similar to steady-state source potentials plus series of terms involving 
various powers of R. By grouping the terms in equation (l4) with 
respect to powers of' I, the following form of the source potential 
convenient for later use is obtained: 


0 N (x,y,z,t) 



aQ2^ + 


+ 


Pm -3 

a 0m R + 


. . .) + l^a^ J + a 12 R 


T,2m-3 

+ . . . + a, R +. 
Im 


, • • . • 


t n/ i 2m-3 \ "It 

s ( a nl r + a n2 R + • • • + ann^ + • • *J + * • J dT l d5 


(17) 


With the terms of the series grouped in this manner, in view of the 
fact that the differential equation (l) is independent of I, it is 
apparent that the coefficient of each power of I in equation ( 17 ) 
is a solution to the differential equation. 

If the values of in equations (11), (12), and (13) are 

put into either equation (14) or equation (17), the Integrations of each 
term can be easily carried out in closed form. Moreover it can readily 
be shown that, when all the terms involving c5 up to a given power are 
taken into account, the differential equation (l) is satisfied to the 
highest power of 55 considered. The boundary condition of tangential 
flow as expressed in equation (4) is satisfied exactly and does not 
depend on the order of 55 considered. 

Putting the values of W(t,T|) in equations (ll), (12), and (13) 
successively into either equation (l4) or equation (17), carrying out 
the indicated integration, and setting z = 0 yields for the velocity 
potential, to the third power of o3 at the upper surface of the wing, 
in the purely supersonic region: 

= - jK^l + Vafg + ^3) (18) 
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where ' 




f l - x - f ^ ^ ( 2 » 2 + 3)* 3 


to ■= x - $ x 2 - -$L (2p 2 + 3 ) x 3 + Ji! (2S 2 + 5)x* 




12M C 


48m 2 


> (19) 


f 3 = | (x - 2 xq ) " (x - 3x 0 ) - -^2 ( 2p2 + 3)x 3 (x - 4x 0 ) 


Derivation of 0^, 0jg, and 0<p^ 

For convenience in the derivation of the velocity potentials in 
the regions of mixed supersonic flow the coordinate system is chosen 
with the origin at one tip and the y-axis coincident with the wing 
leading edge. (See fig. 2(h).) Region T 3 (shown in fig. 1(a)) exists 
whenever the Mach lines from the tips intersect one another on the wing 
(that is, when A(3 < 2). When the Mach lines from the tip intersect 
one another on the wing and in addition intersect the tips ahead of the 
trailing edge (when AP < 1), other regions, not considered herein, have 
to he taken into account, and the determination of the velocity potential 
becomes very cumber some . This discussion is restricted to the condition 
that the aspect ratio he greater than or equal to l/p, that is, AP ==? 1. 

Because of the similarity of conditions in regions Tq and T2, 
the potentials 0Tj_ and 0 r jV3 will he of the same form. It is there- 
fore only necessary to derive one of these, say 0^. The other 

potential can then he obtained by a simple translation of variables. 

The potential 0 t 3 -> as 6i ven later, is a linear combination of 0£j, 0^, 

and 05>2 • 

As pointed out in reference 1, moving doublets may be used to 
satisfy the boundary-value problem in the regions of mixed supersonic 
flow. The potential of the type of doublet required may be obtained by 
partial differentiation of the potential of a moving source with respect 
to the direction normal to the plane of the wing (in the present 
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notation with respect to z). If the expanded form of the unit source 
potential appearing in equation (17) is partially differentiated with 
respect to z, the following expanded form of a unit doublet potential 
is obtained: 

( a 01 I + a 02 R + • • • + a CtaP^ m 3 + * * ') + ? ( a H I + 

_ _ pm -3 \ „n/ 1 

a-^^R + • • • + a -| j^R + • • .)+. . . + s l&ni — + a^/^R + • • • + 

a^^2m~ 3 + . . . . TJdT) d| (20) 

Since the coefficient of each power of I in equation (20) satisfies 
the differential equation and since the differential equation is linear, 
it is permissible to weight these coefficients differently and therefore 
write the velocity potential for region Tq as 

D 0 (^l ff + a 02 R + • • • + a 0n? 2m 3 + • • •) + 


/ 1 2m-3 \ 

+ • • • + a^jj^ +. • . ) 5 + ; • .+ 

R n^ a nl 

where Dq(£,T|) (i = 0,l,2,...n) are the weight factors or distribution 
functions and where the region of integration r is the portion of the 
wing lying in the forecone emanating from the point (x,y,z). (See 
fig. 2(b) for the case z = 0.) 

The problem is now resolved to that of determining the distribu- 
tion function for each series in equation (21) so as to make this 
equation satisfy the boundary conditions. Tne determination of these 
functions can be made with the aid of analogy with known steady-state 
solutions. Examination of equation (21) shows that the first term of 
each series has the form of a steady-state doublet potential and, as 
discussed subsequently, is therefore .a type of singularity that can be 
superposed to satisfy certain conditions of tangential flow for the 
antisymmetrical type of problem in regions of mixed supersonic flow. 

The second term of each series, with the indicated differentiation 


P + + 


+ &nmR 


,2m-3 


.n 


dn dl (21) 




02 = 


ncot 
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performed, is a singular “type of potential, (steady-state source type)}, 
however, it is multiplied by z and would therefore yield zero vertical 
velocity at the xy-plane . All other terms appearing in ea ch series in 
equation (21) are similarly multiplied by z and are nonsingular } 
they therefore yield zero vertical velocity at the xy-plane. In order 
to satisfy the condition of tangential flow it is therefore necessary 
to consider only the first term in each series. 

Retaining, as necessary, only the first term in each series in 
equation (21) and imposing the condition of tangential flow as expressed 
in equation (2) results in an integral equation for the distribution 
functions Dp (i = 0,1,2, .. .n) as follows: 


w(x,y,z,t) 



{^O^l f + D l a ll | + • • • + 


( 22 ) 


Consider the following integral which is of the type appearing in 
a representative first term in equation (21) and which represents, 
as previously mentioned, the potential of a doublet distribution for the 
steady-flow problem: 


lim 
z — >0 


5 icot ^2 






D n a nl 1" + * * • ) dT l 


^ff BK R - dT) (K>0) ( 23 ) 


In steady flow a distribution of this type is convenient for treating 
the problem of satisfying the condition of tangential flow for a dis- 
tribution of normal velocity at the wing surface prescribed independent 
of y^ but proportional to The required integral equation for the ' 

distribution function D^- in the steady case is 

s 2 r r e K 

lim — 2 / / % R d ^ dT l = Cx K (2k) 

Z— >00 oz J J 

where C is a known constant. 

As noted in the following, the distribution function D*- may be 
known provided D 0 is known. The value of D 0 in equation ( 2k) is known 
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for a variety of wing plan forms* In the case of rectangular wings the 
expression for Dq is derived in reference 1 and found to he 


D 


0 


2 

P 


\/Pt)(£ - Pt]) + I sin -1 



(25) 


This expression may he considered as a key for finding any other dis- 
tribution function Dk (K > 0 ) required to satisfy equation (24) for 
the rectangular wing. It follows hy direct substitution and reduction 
that any Dg (K > 0 ) for this wing may he written in terms of Dq, 
equation ( 25 ) ) a s 

% = K! D 0 (d£) K ( 26 ) 

J 0 d 0 d o 


The boundary conditions (a) and (h) can he shown to he satisfied 
hy equation ( 22 ) after the values of Dj£, given in equation ( 26 ), and 
the values of W(£,t]), defined hy equations (ll), ( 12 ), and (13), are 
substituted into equation ( 21 ) for the velocity potential and the 
identity in equation (l 6 ) is utilized. The velocity potential 0p^ 

is thus determined hy equation (21) and these substitutions. At the 
upper surface of the wing the velocity potential 0 ^ is given hy the 

expression 




P 2 cd 2 „ 

M 2 > 2 " M 2 3 


Vd 


2*1 - feiaS + ^ x - ^|)f 2 - ^ F 
' M 2 2M 2 / M 2 


3 + ( 2 P 2 " 1 ) F 4 


6 m 2 


a 


2(x - x 0 )F 1 - 


7Tn2 p 

2 + 2 ico(x - x 0 ) + — p (x - 2 xx 0 ) 
— 2 IT 


Fq + 


P 2 tu 2 f 

2 in) - k__ ( x - Xq ) 


F 3 + °^-(2V 2 + 1)F 4 J 
2M 2 ’ 


(27) 


where the terms are grouped conveniently hy the definition of F n in 
the following integral: “ 


Fn = r ^ 

J o 


sin 


-1 


P ax 

I X 


(n = 1 , 2 , 3 , 4) 


( 28 ) 
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(The functions F n , given in equation (28), and certain related func- 
tions are of particular importance in the remainder of this development. 
Integrated values of this function for the first few values of n and 
expressions for related functions needed later in this, analysis are 
given in the appendix.) 

Examination of equation (27) shows that along the Mach line x = 0y, 
separating region Tq from region N, the expression reduces to 

the expression for given in equation (l8). 

The corresponding potentials for regions T 2 and T^ can now he 
obtained. The potential 0^ is obtained from equation (27) by merely 

substituting 2s - y for y in equation (28) so that 

Gn = xn_1 Gx (n = 1, 2, 3, k) (29) 

Jq 

The potential in region T3 (that is for 1 Aj3 < 2) is a simple 
superposition of the potentials for regions N, Tq, and Tg, as in the 
steady case (see, for example, reference 5 ) , and. may be written as 

0T^ = 0Tq + 0T2 " 0N (30) 


Forces and Moments 

Two types of force and moment coefficients will be derived. First, 
in order to gain some insight into the over -all effect of aspect ratio 
on the forces and moments, expressions for total force and moment 
coefficients are derived. Then, in order to present expressions that 
are more suitable for use in flutter calculations, expressions for 
section force and moment coefficients for any station along the span 
will be derived. 

Total forces and moments . - The local pressure difference between 
the upper and lower surfaces on the wing may be written 



In order to derive expressions for total forces and total moments it is 
only necessary to consider the velocity potential in two regions; either 
regions N and Tq or regions N and Tg. Therefore the expression 
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for the total force, positive downward, on the wing may he written as 
P-= -2 J J Apjj dy dx - 2 J Ap Tl dy dx (32) 


T n 


where Apjj is to he calculated from equation (l8) and the integrations 
in the first term are to he extended over the shaded portion of 
region N shown in figure 2(a), and where App^ is to he calculated 

from equation (27) and the integrations in the second term are to. he 
extended over region Tq. (The integrations in the first term are 
simple and may he performed hy inspection. Those in the second term 
may he readily performed hy making use of the relations given in the 
appendix. ) 

After the indicated integrations have heen performed and all 
position coordinates involved have heen referred to the chord 2h (hut 
the original coordinate symbols maintained), the results can he written 
as 


P = - 8ph^V^k^ Ae lcut 


17 ( L ! + iLg) + cc q (l 3 + i 


(33) 


where the reduced frequency k is related to do and (a hy the 
relations 


i hcD hg^ _ 


m 


and where 

5 ■ *(» ■ H*) 


— 1 M^k 1 

2 ~ pk " p5 " 2A 


3 2 k 33 


^ (4 ♦ P 2 ) 


l 3 = ~2 
3k 2 


^(3 + 3 2 + 63 2 xo) + 


(35a) 

(35k) 


kp 2 x 0 (2 + 3 2 ) 


(35c) 
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L lj. = ~5“(P 2 - 1“- 2P 2 Xq) 

3 J k 


♦ %5 + P 2 

63' 


+ W^Xq) 



I^k 

— g(8 + 43 2 + 203 2 xo + 53^0) 
53 _ 


(354) 


The quanities Lj_ (i - 1,2, 3, 4) are the in-phase and out-of -phase 
components of the total force coefficients, and L 3 being the in- 

phase and L 2 and L^, the out -of -phase components. It will he noted 
that L-]_ and L 2 are associated only with vertical translations of 

the wing and are independent of axis -of -rotation location xq. The 
components L 3 and 14 are associated with angular position anrl 

rotation of the wing about any axis x = xq and depend partly on the 
location of xq. 


The total moment, positive clockwise, on the wing about the 
arbitrary axis of rotation x = x Q is 



- x 0 ) Ap^ dy dx 


~ 2 I ^ “' X °) Ap T 1 dy 


dx ( 36 ) 


If steps similar to those required to obtain equation ( 33 ) are performed, 
there is obtained 


M = -8pb 3 V 2 k 2 Ae IKOt 
a 



iMg) + a Q (M 3 + iM^) 


where 


(37) 


Mi = ~(2 - 3x 0 ) - ^-jr [33 2 + 1 - 4x 0 (2 + p 2 )l 
33 6A3 L U 


m 2 - Jjd - 2xo) - %3 - 4*o> 


235 



r 

1 

1 ( 

3A 

3 2 k 


(2 - 3xq ) - 


^[>(3 2 + 4) (4 - 5xo)] 
53° 


( 38 a) 


( 38 b) 
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m 3 - ^ - ^o) 1 


2p5 


3 + p 2 \+ 4x 0 (p 2 - l) - 8p 2 x 0 2 


JL_ 

3A 


~-p( 2 - 3xo ) 

p 2 k 2 5p b L_ 


^ + 3P 2 ) + 5*o(3P 4 + 3P 2 - 4) 


20p 2 xo 2 (3 2 + 2) 


= 


2 ji(p 2 - 1) - 3Xq(23 2 - 1) + 6p 2 x o 2 ! + ^E_(20 + 4p 2 


3P 3 k 


(38c) 
25 xq + 


403% - 603 2 xo 2 ) + + %(3 2 - l) - 63 2 *o5 - 

gs[20(2 + p 2 ) - 24 xq (2 - 33 2 - p 4 ) - 30P 2 x 0 2 (4 + p 2 []| (38a) 


The quantities 0£ and Mq are , respectively, the in-phase and out -of - 
phase components of total moment coefficients about the axis x = xq 
associated with vertical translations of the wing; M3 and M4 are the 
corresponding components due to angular position and rotation of the 
wing about x = xq. 

It is of interest to note in equations (35)' and (38) that the 
components Lq and Mq do not involve the reduced frequency k. The 

effect of frequency on these two components comes from terms involving 
the frequency to the fourth and higher powers; but for values of k 
thought likely to be encountered in supersonic flutter (k <0.l), the 
contribution of these higher-powered terms to any of the components in 
equations (35) and (38) is, for the most part, negligible. 

Section forces and moments .- The section forces and moments at 
any spanwise station are derived by integrating, the pressure difference 
along the chord for the forces, and the pressure difference multiplied 
by a moment arm for the moments. Since the distribution over the entire 
wing is symmetrical with regard to the midspan section, it is only 
necessary to derive expressions for the forces and moments at any 
station of the half -span adjacent to the origin. (See figs. 1, 2, 
and 3.) 
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Under the restrictions previously stipulated, two cases that can 
arise are considered (see fig. 3): (l) the Mach lines from the tips 

do not intersect on the wing (or A3 > 2) and (2) the Mach lines 
intersect on the wing but the Mach line from one tip does not intersect 
the opposite tip ahead of the trailing edge (or 1 A3 ^=2). Only the 
final forms of the section force and moment equations are given. These 
forms are easily calculated by deriving the pressure difference for the 
different regions from the appropriate velocity potential, makin g use 
of figure 3 to determine the limits of integration for the regions 
involved, and Using the relations given in the appendix to carry out 
the more troublesome integrations. The integrated expression for any 
region can then be reduced to the forms 


P = -4pbV 2 k 2 e iait 


ho 

^r-(Li + iLg) + a 0 (L3 + iLij.) 


and 


= -4pb 2 V 2 k 2 e ljX> ^ 


1*0 

^"(Mp + iMg) 


+ a Q (M3 + iM^) 


(39) 




where the position coordinates are referred to the chord length 2b. 


The components of force and moment coefficients for the half -span 
adjacent to the origin are as follows: 

Case 1 (see fig. 3(a)): For any section between the tip and the 

point where the Mach line intersects the trailing edge, or 

where 0 < y < i. , the components of section force coefficients are 




(^la) 


L 2 “ 



r 


1 

3* 

2k 


•2,. i — 


Mk 
3 


(23 2 - l) Fg - 33 2 f 3 


(4lb) 


L -A 
L 3 - 3 * 


1 - 2(1 - 2x 0 ) W ± + 


2k 2 

2(23 2 + 1)xq = 


r r 


— Fo 


f 2 - 


^3 2 ± 9 i 


(lj-lc) 
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= 


Pit 


1 

"k 

\ 


(2 - xq) Fi - 


W. 3P1+1 a 


f 2 


M 2 k 
+ ^ 


3P | 


3(1 - P 2 + 20 4 + 


2p% - 40%) f 2 - 60 ^ (4 - 3xq) f 3 - (l - T0 2 - 200 


4)f 3 


(4ld) 


where F n (n = 1, 2, 3 , 4), given in the appendix, is obtained from F n , 
equation (28), when x = 2b. For any section between the point where 
the Mach line intersects the trailing edge and the midspan, or 


where 



the components of section forces are: 


_ 

L 1 " 03 

_ 1 M 2 k 

^ " 0k ’ p5 

l 3 = - -^[(3 + P 2 ) + 60% 

J 0k 2 3P 5 L ■ 




> (42) 




60 ?X 


+ 0 £ 


120 2 Xq t 


J 


As a check on the results in equations (4l) and (42) the expressions 
in equations (4l) reduce to those in equations (42) when y = 


The conponents of section moment coefficients for case 1 are as 
follows : For 0 < y < 


% = -Jk 


(1 - 2x 0 ) F x + 2x 0 


20 2 + 1 


P 2 


30 ,t_l ? 




(43a) 




4 

P« 

4p 2 ± i f ; 

p 2 


M 2 k 


m - + s 




(2 P 2-1)(1-^o lf 2 + &1 ^ 3 


P 2 


(43b) 
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M 3 


AM. „ W A ^ - 3xo + 3xo 2 ) = 
k 2^2 " WU Fl + 


(6p^ + 3P 2 - l)(l - 2xp) + ^ 2 xq 2 (i + 2g 2 ) = 


Fg + 


§il ± P 2 ) „ f 

a 2 *0*3 


20P^ + 213 2 + 3 = 


1% = 


k _2 

P*U 


(i - .o) 2 r x + i6i±i ^ - ai±i I 


;2 3 


4l g k f 1 


(] p4 “ 2p2 + 1 )( 1 * - 3x o + 3x 0 2 ) + 1 - Sxo^Jf. 
(2 - 4xo + 3xo 2 ) F3 - ^ ( 83^ + 43 2 - l)f% + 


^(5 p2 + 3)f^ 


and for | < y S 


M X = "^(2 - 3x 0 ) 

3P 3 

M 2 = i(l - 2xo) - g|(3 - 4xo) 




3 = p ^ 2 fl " 2 x o)-- “5C3 + P 2 - 4 x 0 + 43 2 x 0 - 83 2 x 0 2 ) y 

M 4 = ^H Cp2 " 1} ' 3x 0^ 2 P 2 " !) + 632x0 s ] + 

Y f(5 + P 2 ) + 5 xq( 8 p 2 - 5) - 603 2 xo 2] 


3P J k' 

M 2 k 




(^3c) 


(^3d) 


m 



22 


NACA TN 2064 


The expressions in equations (43) reduce to those in equations (44) 
when y = i. The expressions in equations (42) and (44) correspond to 

the more exact two-dimensional components of force and moment coeffi- 
cients derived in reference 4. For values of k< 0.1 these expres- 
sions yield, for the most part, values that are in good agreement with 
those that may "be calculated from the tables in reference 4. 


Case 2 (see fig. 3(b)): For any section between the tip at y = 0 

and the point where the Mach line from the tip at y = 2s intersects 

the trailing edge ^ or where 0< y Is A - the components of section 

force coefficients are given by equations (4l) and the components of 
section moment coefficients, by equations (43). For any section between 
the point where the Mach line from the tip at y 

trailing edge and the midspan ^or where A 

of section force coefficients are 


-i< 


= 2s intersects the 

, the components 


Lq = 


P* 


(Fl + G-l) 


1 +- S A? Z 


P‘ 


+ G2) 


*4 

p 3 


(45a) 


L 2 ' + 


G l) 


M 2 k 


+ |lf[^ p£ 


1)(F 2 + G 2 ) - 3P 2 (F 3 + CJ3O 


( 1_ M g k \ 

\ Pk “ p5 ) 


p5 


(45b) 


t - Jt 
l 3 3* 


1 -(Fq + Gq) - 2(l 


2k c 


2xq)(Fi + Gq) + 


6p^ + 3P 2 - 1 ,= = . ?( g P 2 + 1 ) x 0 ,= = v 

— — (^2 + G2) ^2 (-^2 + G 2 ) - 


6p2 + 5 rf 


(F 3 + G3) 


_ J 

H 

1 

H 

\ 


Pk 2 3P 5 1 

V. 


(3 + P 2 ) + 6 p 2 xqJ 


(45c) 
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L 4 = [ 2 - x o)^l + G l) - ~ - " p 2 1 (F 2 + G 2)] + 

“||3(1 - 3 2 + 2P 4 + 2P 2 x q - 4p lf x 0 )(F 2 + G 2 ) - 


6f3^ - 3x 0 )(F 3 + G3) - Cl - if - 20P 4 )(F4 + G k ) 


-y- |(P 2 " !) " 2 3 2 Xq] + ^|(5 + P 2 + 12P 2 Xq) 
P J k L “ 6p ' ' y 


(45d) 


where . G n (n - 1, 2, 3 j 4) } given in the appendix, is obtained from G n , 
equation { 29 ), when x = 2b. The corresponding components of section 
moment coefficients are 


Mi = 


Pit 


(1 - 2x 0 )(F 1 + G-l) + — ° - ( g ^ * ^ (f 2 + Gg) 


33^1 ( f = 3 + -) + _|_ (2 . 3xo) 


(46a) 




: ji 1^2 + ° 2 ) - xoc^i + + 


M 2 k 


(2p2 - 1)(1 - 2 xn) 

(*2 + g 2 ) + 6xq(F 3 + g 3 ) - 


P 2 


~~t 1 (n + %) 


ik (1 - 2xo) - ~ (3 - 4xo ) 


pk 


2p5 


(46b) 
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(F 2 + G 2 ) - Xo(F x + oj] - 4 ^ (F x + Gjj + 


(6^ + 3P 2 - 1)(1 - 2x 0 ) + 4B 2 Xq 2 (1 + 2B 2 ) _ = 

f (F2 + G 2 ) + 


4 

l 

Bit] 

k2 


6x 0 (l + B 2 ) ,= = N 

(F 3 + G 3 ) - 


1 2 f (f4 + &>; 


- 2i0) - 55 


{(3 + B 2 ) - 4xq(1 - B 2 ) - 8|3 2 xo^] 


(46c) 


M 4 = 

4M 2 k 

„2 


- o 2 

P _ = 2B +1 ,= =• % 2B + 1.— = . 

(1 - xqJ^CFx + Gq) + ~£ — x o(F 2 + G 2 ) - ^2^( F 3 + G 3^ 


1 

6? 


kn4 1_ 


fa-B^ - 2B 2 + 1)(1 - 3 xq + 3xo 2 ) + 1 - Sxo^jCFg + G 2 ) - 


(2 - 4xo + 3xo 2 )(F 3 + G 3 ) - + 4b 2 - IHFIj. + Glj.) + 

6b ^ 


^(5B 2 + 3)(% + g 5 )|J _< ^-|s(b 2 " D - 3x 0 (2B 2 - 1) + 6B 2 x 0 2] + 


3B J k 


M 2 k 


15B 


7 


[ 4(5 + B 2 ) + 5 xq( 8 b 2 - 5) - 60 B 2 x/] 


(46d) 


For the limiting condition of case 2 , that is, when the Mach line 

from one tip intersects the opposite tip at the trailing edge, or A = 

(see fig. 3(c)), the components of section force coefficients are given 
by equations (45) and the corresponding components of moment coeffi- 
cients by equations (46). 


SOME PARTICULAR CALCULATIONS AND DISCUSSIONS 


From the expressions for total force and moment coefficients, 
equations ( 35 ) and ( 38 ), respectively, the over-all effect of aspect 
ratio on the magnitude of the forces and moments can be calculated 
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for particular values of the parameters M, k, xo,,and A. Examinat ion 
of these equations shows that varying some of the parameters might 
cause some terms in the equations to vanish and change sign. For 
example, if xq is continuously increased from some value less than l/2 
to some value greater than l/ 2 , the first terms in the expressions 

for M 2 and M 3 vanish at xq = ^ and change sign when xo "becomes 

greater than l/2. In particular, decreasing the aspect ratio decreases 
the components of force and moment coefficients Lp, T^, X 3 , Mq, M 2 , 
and Mj "but increases the two important components L]^ and M^. 


Although the effect of aspect ratio may change considerably with 
only a small change in any one (or more) of the parameters M, k, 
and xq,. some insight into what the over -all effect might be can be 
gained from calculations of all the components of total force and 
moment coefficients for various values of M and A and fixed values 
of the parameters k and xq. Results of such a set of calculations 
are presented in figures 4 to 7 - 


In figures 4 and 5 the components of total force and moment 
coefficients for various values of A and for xo = 0.4 and k = 0.02 
are plotted as functions of M^. The curves in these figures calculated 
for infinite aspect ratio correspond to the two-dimensional results of 
reference 4. The dashed curves represent calculations for aspect ratio 
and Mach number combinations that cause the Mach lines from one tip to 
intersect the opposite tip at the trailing edge so that along the 
dashed curves the aspect ratio is not constant but varies with M^ 
according to the previously given expression 



The difference, at any value of M^, between the dashed curves and the 
curves corresponding to infinite aspect ratio in figures 4 and 5 is, 
therefore, for the chosen values of k and xq, the maxi mum effect of 
aspect ratio on the components of total force and moment coefficients 
for a rectangular wing under the restrictions of the foregoing analysis. 
It will be noted in figures 4 and 5 that, when the aspect ratio is 
small, the deviation of the three-dimensional results from two- 
dimensional results may be quite large. 

In figures 6 and 7 the components of the total force and moment 
coefficients are plotted as functions of aspect ratio for Xq = 0.4, 
k = 0.02, and some particular values of M. It will be noted in these 
figures that all the components of force and moment coefficients undergo 
rapid changes with respect to varying aspect ratio when A becomes less 
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than 4 or 5- It may "be remarked that the directions of the changes 
■with respect to aspect ratio appear to he such that they would have 
favorable effects on the flutter characteristics of a wing. 


The spanwise distribution of the components of section force and 
moment coefficients computed from equations (4l) to (44) for A = 4, 
xo = 0.4, k = 0.02, and M = 2 are plotted in'figures 8 and 9. The 
portions of the curves in these, figures corresponding to value's of y 


in the range i < y ^ A - i are the two-dimensional values and the 
effect of aspect ratio may be noted in the tip regions, 0 y = i 
and A - = y ^ A, as deviations from these two-dimensional values. 


In conclusion, it may be stated that, in regard to the effect of 
aspect ratio on supersonic flutter, an important item that has not been 
discussed herein but can be studied for any particular case with the 
aid of equations (33) and ( 38 ) is the change in center of pressure, 
associated with prescribed motions of the wing, with change in aspect 
ratio. An investigation to find the effect that thickness might have 
on the center -of -pressure location is also needed. An extension of the 
foregoing analysis to include the effect of an aileron as an additional 
degree of freedom would follow in a straightforward manner. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va. , January 6, 1950 
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APPENDIX 

SOME INTEGRATED VALUES OF F n , Gn, f n , AND 
OTHER RELATED FUNCTIONS 


Values of F n and Gn . - The values of the functions F n , equa- 

^-*- on (28 ) , and Gn, equation (29), for the first few values of n are 
as follows : 

nx 

F n =J o x n_1 sin -1 \fiy/x dx (n = 1, 2, 3, . . . ) 

F 1 = ^(x - Py) + x sin -1 \/3y/x 

F 2 = ^'" + ^ 2f3y \/Py( x “ 3y) + \ sin -1 \J&yjx 

f 3 - ^ ^ vSTTT^) + $ sln -i ^ 

% - g * 3 ' ^ * l6ey fib- Br) ♦ ^ »*»■* 

t? _ 35x 1< ' + 4opyx3 + 48(3 2 y 2 x 2 + 64p3y3 x + i28sV^ IZ~, : 

5 1575 V Py(x " Py) + 

5 

-5- sin' 1 VPy/x 


G n = f Q xU ~ 1 s 111 " 1 \J ~ 2S X ~ y) ^ (n = l, 2 , 3, 

Gi =^P(2s - y)£ - p(2s - y)J + x sin" 1 ^ 25 ^ 2? 


• •) 
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xF n dx - — (x F n , - F n+ 2 ) 

x 2 F n ax = |(x 3 F n + F n+3 ) 

x 2 - 2 xxq i 

(x - Xo)F n dx = g F n + x O F n+l “ 2 F n+2 

O x3 - 3x 0 x2 + 3 x 0 2x 

(x - X 0 ) 2 F n ax = 3 F n - Xq F n+1 + xoF n+2 


|^n + 3 


Corre sponding integral relations for G n may be obtained, from these 
relations by simply replacing F by G. 

Integral relations for F n . - Integral relations for F n that may 
be used, in calculating total forces and moments from sectional forces 
and moments are as follows: 

- l x dy » 

,< /s * * - 5 

2h / 0 1/P f 3 ^ “ if 

r !/P = 3brt 

Zb J 0 F 4 = h 

2b J 0 F 5 ^ “ fefT 
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(d) L k . 

Figure 4.- Concluded. 
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Figure 5.- Continued. 






Figure 5.- Continued 






Figure 6.- Components of total force coefficients as functions of A for Xq =0.4, k = 0.02, and 

various values of M. ' 9 
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Figure 7*- Components of total moment coefficients as functions of A for Xq = 0.4, k = 0.02 and 

various values of M. ’ 
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Figure 7 .- Continued 
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